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SUMMARY 


The objective of this study was to theoretically examine the open- 
loop combustion response of monepropellant droplets and sprays to im- 
posed pressure oscilxationso The theoretical model is an extension of 
an earlier response model developed for one— dimensional liquid monopro- 
pellant strand combustion^ The theoretical approach employed pertur- 
bation analysis through first order, yielding linear response results » 

The major assumptions of the analysis were; a one-step, irrevers- 
ible gas phase reaction; phase equilibrium at the liquid surface; con- 
stant liquid, phase physical properties; variable gas phase physical 
properties; and negligible forced convection^ The latter assumption 
is known to be valid at high reactivities, where the reaction zone is 
close to the liquid surfaeec Transient liquid phase analysis was 
coupled with both transient and quasl-^steady gas phase models. The 
calculations emphasized the quasi-steady gas phase eassj since un- 
steady gas phase effects largely became important at frequencies 
greater than the most probable frequency range for combustion insta- 
bility, The calculations were conducted using properties correspond- 
ing to hydrazine decomposition. 

The results of the study may be summarized as follows: 

STEABY GQMBUSTION 

lo The zero-order or ste :.dy state calculations gave predictions 
in good agreement with earlier measurements of liquid surface tempera- 
tures and droplet burning rates in combustion gases, for hydrazine. 

2. Small droplets essentially evaporate with little reaction 
effects while large drops (about two orders of magnitude larger) burn 
with surface mass fluxes equivalent to liquid strands. A simple ap- 
proximate formula was developed which estimates burning rates for all 
sizes with maximum errors on the order of 10%. 

DROPLET RESPONSE 

1. Small drops give low or negative response at all frequencies. 

2. Large drops provide a response near unity at low frequencies. 
Response declines at frequencies somewhat in^ excess of the character- 
istic frequency of the liquid phase thermal wave. - If temperature 
gradients of sufficient magnitude are present in the liquid phase, 
this deeiine is preceeded by a response peak in excess of unity. 

3. If a response peak occurs at liquid phase thermal wave fre- 
quencies, it is associated with a frequency regime where the surface 
temperature oscillation leads the pressure oscillation. At higher 
frequencies the temperature oscillation lags the pressure oscillation, 
and this behavior is associated with declining response. 



4. Large drops have a second response peak of order unity at 
frequencies associated with the gas phase thermal wave; however, the 
frequency range of this peak is generally higher than frequencies 
associated with combustion instability. 

SPRAY RESPONSE 

1. The total response of a monodisperse injected spray is less 
than the response of the initial droplet size; although the frequency 
range of response peaking is similar. 

2. For a polydlsperse injected spray, the response depends 
strongly on the spray distribution. Even a low percentage of large 
droplets can yield an appreciable response, due to their long life- 
time and high response. 

3. At high pressures and low frequencies in the range associated 
with combustion instability, droplet lifetimes approach the cycle 
period for drop sizes of technological importance. The response 
associated with the lifetime was not considered in this investigation, 
and desenjes further study. 
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CHAPTER I 


INTRODUCTION 


General Statement of the Problem 

The phenomenon of rocket combustioii Ihstability is the direct 
result of the coupling between the combustion and fluid dynamics of 
the overall combustion system. This coupling enables oscillatory 
energy from the combustion process to sustain these instabilities. 

The oscxllatxons will decay only by decreasing the coupling mechanism 
or by Insuring the damping process, either natural or man-made, has 
sufficient magnitude to dissipate the supplied oscillatory energy at 
a faster rate than it is produced. 

Almost every rocket engine development program has encountered 
combustion instabilities in one. form or other; thus, understanding this 
phenomena is of major importance. These oscillations, which have 
frequencies in the range of 500-30000 Hz [1], are capable of 
producing vibrations, increased heat transfer rates, decreased 
performance, uncontrolled impulses, and variations in the thrust 
vectors of rocket combustion chambers [1]. Since rockets often 
carry humans and/or expensive and very sensitive? equipment, the 

understanding and prevention of these destructive oscillations are of 
utmost concern, 

Instabxlities can be initiated spontaneously from the inherent 
noise in the combustion and mechanical processes, or they can be 
actxvated by a finite disturbance such as a chamber pressure 
fluctuation or an explosion of accumulated propellants [1]. 


z 

Combustion instabilities can be suppressed by the use of 
acoustic liners, Helmholtz reasonators , baffles, other damping 
devices, or by the redesign of the combustion chamber, in order to 
eliminate resonant modes. The only true test of the stability of a 
rocket engine, however, is from an actual test firing. This further 
underscores the importance of understanding the coupling mechanism 
between the combustion process and the fluid dynamics. 

The objective of the present study is to investigate one 
aspect of the combustion instability problem, namely, the response of a 
burning liquid monopropellant to imposed pressure oscillations. The 
response of a Single droplet is considered first, followed by an 
estimation of the response of an entire spray. 

Hydrazine is selected as the raonopropellaut fuel in this study, 
since the hydrazine family of fuels are extensively used in today's 
space rocket programs. For example, the Apollo command module 
thruster and the lunar ascent and descent engines use hydrazine 
fuels. The major emphasis of this study concerns the theoretical 
analysis of hydrazine droplet combustion and burning rate response to 
fluctuations in pressure. 

1 . 2 Previous Related Studies 

Consideration of previous work is divided into three categories; 
liquid monopropeilant droplet combustion, quasi-steady monopropellant 
spray combustion, and oscillatory combustion, 

1.2,1 Monopropellant Droplet Combustiori . A detailed summary of 
experimental hydrazine droplet and strand (large drop limit) combustion 
has been presented by Allison [2, 3], 
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Fendell [4] obtained a closed form analytical solution for the 
burning of monopropellants in an inert atmosphere at the zero 
activation energy limit. 

[5] restricted his approach to adiabatic burning with 
a low reaction rate and a one-step overall reaction. He developed an 
analytical solution for small drops having low chemical reaction 
rates and large activation energies. His analysis is limited by the 
adiabatic and low reaction rate assumptions . 

Spalding and Jain [6], Jain [7], Tarifa [8], and Allison and 
Faeth [9] have investigated monopropellant droplet combustion using 
the thin flame approximations . Spalding and Jain [6] and Jain [7]: 
used the thin flame approximation in order to obtain an analytical 
solution for monopropellant droplet burning rates. The energy 
release of the decomposition reaction is assumed to occur in a thin 
shell at some distance from the drop. The reaction rate dependence 
on temperature is quite strong for large activation energies, and 
the flame surface is located where the unreaeted gas flows into the 
thin shell at the laminar burning velocity of the mixture, Spalding 
and Jain [6] were able to estimate the droplet burning rate as a 
function of the laminar flame speed, which is known from experiment, 
for adiabatic combustion. 

Tarifa [8] assumed an Arrhenius global reaction rate in a 
flame zone of finrte thickness, and obtained an analytical solution 
for monopropellant combustion in an atmosphere of inert gases. Allison 
and Faeth [9] used Spadling's [6] thin flame approximation in order to 
relate flame position to the laminar burning velocity. An Arrhenius 
correlation with two adjustable parameters was used to approximate the 
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laminar burning velocity„ Allison and Faeth [9] also obtained 
experimental droplet burning rate results which compare favorably 
with results obtained in the present study. 

Lorell and Wise [10] considered the ease of a single, 
unimolecular overall reaction with no backward reaction^ A two point 
boundary value problem was numerically solved in order to find the 
burning rate eigenvalue. 

Later work by Allison and Faeth [11] employed one— step 
kinetics to describe the combustion of a liquid hydrazine strand. 

The steady state conservation equations were integrated numerically, 
using an asymptotic expansion to formally satisfy the boundary 
conditions at large distances from the liquid surface. This procedure 
does not involve any approximate mathematical treatment of the 
combustion process, and since it minimizes uncertainties in the model, 
this approach was employed in the present study, 

1‘2.2 Quasi-Steady Monopropellant Spray Combustion , The 
problem of combustion instability is very complicated and contains 
many variables. There are several analytical methods to describe 
combustion instabilities in combustion chambers; however, they are 
severely limited to idealized problems, and almost none have been 
applied, to the monopropellant case, A review of the current literature 
on this subject can be found 5.n Reference [1], 

The present theoretical study uses a one dimensional, spray 
combustion model. This model considers several drop size groups, in 
order to simulate actual spray characteristics. The response 
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contribution from each of these groups is summed to obtain the overall 
spray response, 

1.2.3 Oscillatory GombustioU c A description of much of the 
available literature on combustion instability is given in Reference 
[1], The present summary primarily considers liquid burning rate 
response to imposed pressure oscillations. 

Crocco [12, 13, 14] uses the sensitive time lag theory in which 
a time lag exists between the time a fuel particle is injected into 
the combustion chamber and the time it is burned to form the final 
combustion products. With this approach, the combustion-gas dynamic 
coupling process is not explicitly considered, and the time lag 
can only be found from experimental observations, 

Priera [15, 16, 17] realized that the thermal energy and the gas 
release rates are of utmost importance in analyzing the combustion 
process, and he numerically solved the gas dynamic equations for the 
combustion chamber to determine these factors. Fluctuations in the 
surrounding gaseous flow field and the unsteady burning rate response 
of the droplets, which is the source term in the gas dynamic equations, 
were also investigated by Priem, 

Williams [18] and Strahle [19, 20] used a bipropellant fuel 
and assumed the Eurke-Schumann thin diffusion flame model in the 
leading edge boundary layer of the burning droplet. The droplet was 
subjected to longitudinal standing wave fluctuations while the 
liquid temperatures were assumed to be constant. The response curves, 
which depend on frequency, were flat with a peak at extremely high 
frequencies. The difficulty with this analysis arises from the fact 
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that liquid temperatures are not constant for oscillatory combustion, 
even at low frequencies, while the gas phase can be considered to be 
quasi-steady except at very high frequencies and pressures [21]. 

Heidmann and Wieber [22] also investigated bipropellant 
droplets, assuming that at each point in time, the gas phase 
combustion rate equalled the vaporization rate of the liquid fuelc 
This is not a good assumption at high frequencies where gas phase 
transient effects become important. However, they did find response 
peaks for quasi-steady drops . 

In the bi.propellant analysis of T'ien and Sirignano [23], the 
gas phase was considered to be quasi-steady with respect to the 
liquid phase at low frequencies. The peaks in the response function 
were caused by the liquid phase thermal lag; however, the response 
was not large enough to provide a mechanism for instability. While 
the liquid phase thermal wave does not give sufficient peaks for 
bipropellants , this may not be Che case for monopropellants , which 
were not considered in Reference [23] . 

Allison and Faeth [11] used a monopropellant hydrazine strand 
model, similar to a model used by T'ien [24] for a solid propellant 
response study. They assumed low pressure phase equilibrium atthe 
liquid surface, unsteady liquid phase, and both quasi— steady and 
unsteady gas phase. The liquid phase thermal wave for strand 
combustion yielded peaks capable of causing instability. Liquid fuel 
response was measured at frequencies representative of the response 
of the liquid phase thermal wave. Agreement between theory and 
experiment regarding both the magnitude and phase lag of the combustion 
response was found to be good. The next logical step in this analysis 


7 

would be the investigation of oscillatory combustion of droplets. 

The present study seeks to complete this extension. 

1.3 Specific Statement of the Problem 

The development of a mathematical model of oscillatory 
combustion of liquid monopropellant droplets is of significant value 
due to the application of these propellants in actual rocket engines. 

Very few, if any, experimental studies of nonsteady liquid 
droplet combustion have been conducted. The past procedure has been 
to obtain pressure-time data from actual rocket engine tests and then 
attempt to deduce burning rate response. 

Allison and Faeth [11] studied steady and unsteady burning 
of liquid strands of hydrazine, and were able to obtain very good 
agreement between experiment and theory. These one-dimensional 
strands approximate the large drop limit, since the curvature effect 
of a large spherical drop is negligible. 

The present study will extend the work of Allison and Faeth 
[9, 11] to derive a mathematical model for unsteady oscillatory droplet 
burning. The large drop values will be compared with the one- 
dimensional strand burning results of Allison and Faeth [11]. The 
forcing function will he the gas phase pressure oscillations, and 
the burning rate or open-loop response will be examined. 

Thus, the specific objectives of the present study are; 

1. Develop a mathematical model for the burning of a liquid 

monopropeliant droplet for both steady state and oscillatory 
pressure conditions. 
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2, Theoretically investigate the steady state droplet burning 
for both wet and non-wet bulb conditions, and compare the 
results with available measurements « 

3c Investigate the theoretical burning rate response for 
wet and non-wet bulb conditions c 
4c Develop a mathematical model for one-dimensional, quasi- 
steady spray combustion, and calculate the overal spray 
combustion response for both wet and non-wet bulb 
conditions. 



CHAPTER II 


THEORETICAL CONSIDERATIONS 

2,1 General Model and Assumptions 

The theoretical model under consideration involves a droplet 
that is burning in the a.bsence of convection, yielding a spherically 
symmetric flow field- I^hile convection is important under actual 
combustion chamber conditions, there is evidence to indicate that 
the region where the droplet has a low (or zero) relative velocity 
with respect to the aiibient gas , is a very critical zone in 
determining combustion instability characteristics [1], For simple 
evaporation without decomposition, the response due to velocity 
effects has been found to be quite low [23], In addition, the 
decomposition process near the droplet surface reduces the influence 
of convection over a fairly wide range of conditions [25], [26], 

This occurs since convection only affects the process when the outer 
edge of the flow field interacts with the reaction zone, a situation 
that is limited to very weak reactivity or very high Reynolds 
numbers. Therefore, neglecting convection for monopropellant droplets 
puts f e^jer limitations on the practicality of the calculations , than 
would be the case for bipropellants. 

Similar to earlier response studies [11], [23], the ambient 
pressure is assumed to be oscillating with a wavelength that is long 
in comparison to the dimensions of the combustion field of the 
droplet. The period of oscillation, however, is assumed to be short 
in comparison to the total droplet lifetime so that consideration does 
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not have to be given to large changes in the position of the surface 
during a pressure oscillation,. This assumption allows the analysis 
to proceed while only examining oscillatory solutions. The frequency 
regime where this assumption breaks down will be discussed later. 

Since the combustion rate of monopropellant droplets is much 
higher than bipropellants, the usual assumption of a constant liquid 
phase temperature (in the zero order) is less valid than for 
bipropellants. Therefore, the presence of mean liquid phase 
temperature gradients must be considered in the analysis. Also, 
examination of constant mean temperatures will be ma.de, in order 
to include the conventional steady combustion model, by equating the 
bulk liquid temperature to the wet bulb temperature at a given 
pressure , 

A sketch of the theorecical model is shown in Figure 1, The 
process is examined at an instant of time when the droplet radius is 
r^*. Formally, this radius is taken to be fixed so that the mass 
flux of fuel is time varying. This actually corresponds to porous 
sphere combustion; however, the two cases are equivalent a.s long as 
the density in the liquid phase is large compared to that of the gas 
phase; and the period of oscillation is small in comparison to the 
lifetime of the droplet [27], IJhen this is true, the motion of the 
surface with respect to the mean surface position is negligible and 
can be ignored. Exceptions to this assumption arise near the critical 
point and the present analysis is not valid in this regime. The 
response portion of the analysis is invalid at frequencies having an 
oscillation period comparable to the lifetime of the droplet. 
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In the true model of monopropellant droplet combustion, the 
radius is not constant, but decreases with time, while m^* is zero. 
Also, as the droplet becomes smaller due to gasification, the 
temperature throughout the droplet tends to approach the surface or 
wet bulb temperature. 

For generality, gas phase transient effects are included in 
the model. The effect of variable properties is also included, so 
that the model is equivalent to the earlier strand combustion analysis 
[ 11 ]. 

The remaining assumptions of the analysis are similar to those 
of Reference [11]. They are as follows: 

1. The process is spherically symmetric with a Mach number 
less than unity and negligible body forces. Inertial and 
viscous terms in the momentum equation are neglected. 

2. The reaction process is premixed and laminar, A one-step, 
irreversible chemical reaction takes place in the gas 
phase and any time lags associated with the chemical 
reaction itself are negligible, i»e,, the chemical reaction 
is locally quasi'^steady and obeys an Arrheviius equation, 
Ghemical reaction is neglected in the liquid phase, 

3. Radiation heat transfer is neglectede 

4. The gas phase is taken to be an ideal gas and the Lewis 
number is assumed to be unity , 

5. All gas diffusion coefficients are equal, all molecular 
weights are equal, all gas phase specific heats are equal 
and constant, the gas phase thermal conductivity is 
independent of composition and varies linearly with 
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'-emperature , and tlie liquid is composed of a single 
chemical species having constant properties. 

6, The combustion products are assumed to be insoluble in 
the liquid phase and the gas phase fuel mass fraction at 
the liquid surface is given by the Clausius-Glapeyron 
equation 0 As in the ease of the gas phase reaction, the 
equilibrium at the surface is assumed to occur rapidly 
in comparison to other transient effects in the system. 

7. The wavelength of any periodic pressure disturbance is 
assumed to be long compared with the diameter of the zone 
involving active combustion. Consideration of the momentum 
equation, along with Assumption 1, then Implies that 
pressure is only a function of time. 

A discussion of the applicability of these assumptions is 
provided in References [3] and [11]. 


2 , 2 Gevernihg Equations 

2.2.1 Dimensional Equations . Considering the above assumptions, 
the dimensional conservation equations can be obtained in spherical 
coordinates from the general equations given in Reference [28]= 

Liquid Phase o ^ r* ^ r^* 

Gonservration of Mass 

0 * V * = m.* - constant (3-1) 

fr f 


Conservation of Energy 


2 9T* . J.2 j, 3T* 

o * Q * + p * c r* ■ V * 

*^f ' f 9t* “^f £ 


r 9r* 


I dr 


9T* 

9r* 


( 2 - 2 ) 



14 


where the mass flow rate per unit solid angle, is only a function 

of time. 


Gas Phase r * < r* < “ 
— s — 

Conservation of Mass 




(2-3) 


Conservation of Species 


9^. 3Y. I 9Y, 

P* arl" + P*v * TT-r = -sAr I p*D* ^ 

9t* r 8r* ^,^2 9r* | 3r* 




Conservation of Energy 


i=l, N 


(2-4) 


pAG * ^ + D*V AC A S - -i- ^ 

^ S 3t* ^ \ ^P 3rA ■ ^^2 8 ^ 


© - P - C2-5) 


where the stoichiometry of the one-step gas phase reaction is assumed 
to have the form 


N N 

i p.' [i] £v " [i] 

1=1 1-1 

The other terms in Equations (2-4) and (2-5) are defined as 
Wp* = B*TA 


p*Y ' 

n 

-E*' 

F 

exp 

[MARATA 

RAX* 


V. = 
1 


V. " - V. ’ 
1 1 


F 


N 


= - Z "^1 

i'Al ^ ^ 


where = 
F 


^1. Since the gas is ideal 


( 2 - 6 ) 


(2-7) 

(2-8) 

(2-9) 


pA = p*rat*/m* 


(2-10) 
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The gas phase mass flow rate per unit solid angle is defined 


as 


m* = 


p*r* 


2 

V 


* 


r 


( 2 - 11 ) 


The thermal conductivity is taken to be proportional to 
temperature 


A* - A. A) 

oo ' ' 00 ^ 


The unity Lewis number assumption provides 


Le = 


A* 


p*D*Cp* 


= 1 


( 2 - 12 ) 


(2-13) 


Boundary Gonditions 

Since the liquid is Incompressible, it has a constant 
temperature at the origin or center of the drop , 

r* = 0 T* = (2-14) 

In the case of steady combustion at the wet bulb temperature of the 
droplet, T^* would be equal to the wet bulb temperature^ However, 
since Tnonopropellants have higher gasification rates than most 
hipropellants , cases exist where temperature gradients are present 
in the liquid phase throughout the lifetime of the droplet „ 

Therefore, T^* is taken as a parameter for the analysis in order to 
investigate the effect of these gradients on the response. 

Temperature and pressure are continuous across the liquid-gas 
interface. 


T * 
s- 


^"■Eii;pE,oDUGiBn.rrY of 

OEIGINAL PAGE IS 


2-15) 
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P * = P * 
s- s+ 


(2-16) 


Conservation of mass at the interface reduees to 


P* V * = D*V * 
fs s 


(2-17) 


The conservation of energy eciuation, where the heat flux 
reaching the surface goes into the heat flux in the liquid phase and 
the heat of vaporization of the gasifying fuel, reduces to 

8T*1 


A.* = X* 




3T* 


9r* 


- P/ V *[(G *)T * + L*] 

s+ fs t s 


(2-18) 


The insolubility of the product species in the liquid phase 
implies that the total mass flux of the non-fuel species is zero at 
the surface of the droplets The conservation of species equation 
then reduees to 

9Y, 


p*B* 


9r* 


= Pf*v *(Y^ - 1) 

s+ ^ ""fs 


(2-19) 


The fuel mass fraction is related to the temperature at the 
surrace by the Glausius-Clapeyron equation. 


'^Fs ~ P* 


r -L vv 
f 


s 


( 2 - 20 ) 


Far from the drop surface all the fuel is consumed and the 
temperature is only a function of time . 


rVV ss CO 




T*=T*(t*) 


( 2 - 21 ) 


where T«(t*) is a known function depending upon the specified 


pressure variation. 



17 

2.2.2 Nondlmensional Equations . The previous equations are 
simplified by transformation into the following dimensionless 
variables . 


a = a*/p^* 


r = r*/r * 


D = D*/B * 

' OOQ 


t = t*X */p *C *r 
°°o '^“o p s 


E = E*/R*T * 
oop 


■p = T*/x * 

OOQ 


L = L*/G *T * 

P OOQ 


V = V *0 *G *r *A * 
r r “o p s “o 


= L.*/R*T * 

f f o°G 


6 = G */G * 

f p 


m = m*C */r *A * 
p s “O 


5- = a */a * 
f f “O 


P = P*/Po* 


X = X*/X * 

OOQ 


0 = q*/G *T * 

^ ^ ' P OOQ 


P == P*/P„Q* 


(2-22) 


The conservation equations can now be retJritten in dimensionless 
form as follows : 


Liquid Phase 


r < 1 


f fr f 


f 9 


^ ^ 

9t ^f r 9r ’ 2 Sft 

r 


2 9T 


(2-23) 

(2-24) 


Gas Phase 


r > 1 


3p . 1 3 , 2 . 

17 V = ° 

r 


(2-25) 
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2 

m = pr V 

r 


(2-26) 


3Y . 


+ V 

r 



1 S 




+ V. w 
1 


i = 1, N (2-27) 


1 


3T 

8t 


+ V 


111 

r 3rJ 


1 

2 

r 



X 


3r 



dp 

dt 


+ q w 


(2-28) 


wTiere the dimensionless reaction rate is defined as 


w = A exp (-E/T) 


(2-29) 


and A is a reaction parameter which is related to the drop size 

Q r P * 

A = -£ ^ I ° 


OOQ 


cog 


(2-30) 


Th:e assumptions of an ideal gas, unity Lewis number, and thermal 
conductivity proportional to temperature, reduce to the following 
relationship 


T = X = p0 = p/p (2-31) 


Boundary Conditions 

The dimensionless boundary condition equations are: 
At the center of the drop, r = o 


T = T, 


(2-32) 


At the drop surface, r = 1 


1 = T 

s- s+ 


(2-33) 


P = P , 
s- s+ 

hH| -PV^U1-0)T + L] 

»S- IS+ 


(2-34) 

(2-35) 
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's+ 

= a exp (-L^/T) 

Far from the drop surface r = «> 


T = T(t) 


(2-36) 

(2-37) 


(2-38) 


(2-39) 


The above equations can be further simplified by introducing 
the Shvab-Zeldovieh variable, 0 , a new dependent variable based on 
the fuel mass fraction 


9 = q Yp + T 


(2-40) 


After substituting Equation (2-40) into Equations (2-27) 
to (2-29) and (2-36) to (2-38), along with some minor simplifications, 
the pertinent equations can be summarized as follows: 


Liquid Phase 


r < 1 


Gas Phase 


9T ^ 9T 
8t ^fr 3r 


r > 1 


A 

r2 9r" 9r^ 


9p , 1 9,2.' 

"T 5? ® 


p|i+v (A|i 

dt r dr ^2 dr 9r 


9T 


3.T _ 1 9 / 2„ 3T 


I = ^ I r-T ^ 




dt 


9.t ■ "r 9r ]' 2 9? ] ^ ^ Wj 7 ? + ^ 

w = A q^" p’^ T°^’^(0 - T)“ exp(-E/T) 


(2-41) 


(2-42) 

(2-43) 

(2-44) 

(2-45) 


Rfi^>iiGDiiciBiLrry of the 
origin^al page is poor 
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Boundary Conditions 
r = o 


r = 1 



T = T_ 

(2-47) 


a, 9ll 


1 

37 - 

^17 - P V [(l-g)T + L] 

(2-48) 

/s- 

/s+ ^ 

T 9(6 - 

T) V 


^ 3r 

/s+ " ^ T - q) 

(2-49) 

0 - 

• T = aq exp(-L^/T) 

(2-50) 


0 = T = T(t) 

(2-51) 


Analysis . The imposed pressure oscillations 
can be decomposed into a constant mean pressure and an oscillatory 
component having a small amplitude, e, Transforming the dependent 
variables in the same manner, substituting into the conservation and 
boundary equations, and separating into like powers of e. results in 
zero-order or steady state equations (e°) and first-order or oscillatory 
equations (e^) , The higher powers of e are disregarded in the present 
analysis, which considers only linear response. The assumed 

pressure oscillation and the dependent variables are represented as 
follows ; 
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p(c) = 1 + E 

T(r, t) = T^(r) +ET^(r)e^*^^ 

0(r, t) = 6^(r) +E0^(r)e''-‘*^‘' 

rafr, t) = m^(r) + s 

m^(t) = ifi^^ + E i; 2 - 52 ) 

where T^, 0^, and are complex quantities. 


2.3 Zero-Order Problem (Steady State) 

2„3.1 Equations . Substituting Equations (2-52) into 
Equations (2-41) to (2-51), yields the following zero-order, steady 
state fequations : 


Liquid Phase 


r < 1 


d ( 2 *^^0 ™o _ „ 

dr P dr ^o| - ° 


- Tn = constant 
fo o 


(2-53) 

(2-54) 


Gas Phase 


r > 1 


I 

dr ■ o dr 


- m 0 I = 0 

O Of 


(2-55) 


dT \ 

dr'"^ dJ^ - % ^o r ■ ^o^"" exp(-E/T^) 


= 0 


(2-56) 


Boundary Conditions 


r = o 


T = T. 
o 1 


(2-57) 




frr—TVT 
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r = 1 


X 


f 



T 

o 


s+ 


- A [(i-e)T + L] 
o o 


T 

Q 


d(0 - T ) 

o o 

dr 


s+ 


= III (6 
o o 



(2-58) 


(2-59) 


9 

o 


- = a q exp(-L^/T^) 


(2-60) 


r 


00 


0 = T = 1 

O O 


(2-61) 


These zero-order equations are sufficiently general to allow 
for any ambient conditions, either adiabatic or non-adiabatic« Since 
most combustion chambers are essentially adiabatic, and monopropellants 
do not involve a local enthalpy defect in any stream, a condition of 
adiabatic combustion can be applied to simplify the calculations. (This 


is Only strictly true under the conditions of the present model which 
allows for temperature gradients in the liquid phase. If there is 
bulk, heating of the liquid, small enthalpy defects will arise in the 
spray,) This leads to the zero-order compatability condition between 
dimensionless properties as follows: 

q = 1 - 3 T^ + L (2-62) 


The flame temperature can be computed from the dimensional 
form of Equation (2-62), 


2»3.2 Liquid Phase Solution , An analytical solution can be 


obtained for the zero-order liquid phase Equation (2-53) from the 
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boundary conditions T = T. at r = o and T = T at r = 1. as 

o i o os ’ 

follows: 


T - 

o i 

T - T, 
os i 


L 

) ° 1 

1 , 


r 'W] 


(2-63) 


The boundary condition of Equation (2-58) can now be simplified 
by substituting the derivative of Equation (2-63) at r = 1, and 
noting that B = 


dT 


o 


j , = rii (T - B T. + L) 

o dr / , o o i ' 

's+ 


2,3.3 Gas Phase Solution. Por the case of adiabatic 


(2-64) 


combustion, integration of Equation (2-55) once; application of the 
boundary conditions at the liquid surface; and simplifying through the 
use of Equation (2-62) , yields 

d0 


2 

r T 


^ — m 0 = - m 

o dr o o o 


(2-65) 


Integrating Equation (2-65) again, after transformation to eliminate 
explicit dependent on T^, yields after application of the outer 
boundary condition 



( 2 - 66 ) 


as a consequence of adiabatic combustion. Substituting Equation 
(2-66) into Equations (2-59) and (2-60) yields 



+ T 

o 



(2-67) 


T = 1 - a q exp (-L _/T ) 

os ^ f os 


( 2 - 68 ) 
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Equation (2-68) may be solved with an iteration procedure for the 
temperature at the liquid surface. The outer boundary condition, 
Equation (2-61) is unchanged. 

The zero order problem is now reduced to the solution of a 
single, nonlinear, second order, ordinary differential equation. 
Equation (2-56). Equations (2-61), (2-67) and (2-68) provide three 

boundary conditions to define the problem and determine the eigenvalue, 

ih ^ 
o 

In the general case. Equation (2-56) is nonlinear and must be 
xntegrated numerically. Since any practical range of numerical 
integration must be finite with respect to r, the outer boundary 
condition. Equation (2-61) was handled by considering the asymptotic 
behavior of Equation (2-56) for large r. 

To find the form of Equation (2-56) for large r, let 

(2-69) 


Where ? is a small parameter. The radial coordinate is also stretched 
by defining 


r = £ n 


(2-70) 


where a' is a constant to be determined by proper matching. Keeping 
the lowest order terms, in T^ the equation becomes: 


dn ^ dn 


( 1 ) 


- m £ 
o 


dT 


( 1 ) 


dr) 


j.2a' + n - 1 2 , . ^ (lin 

" ? n A exp(-E) T ^ ^ = 0 

o 

(2-71) 


The form of the asymptotic solution depends upon the reaction 
order, n. The combustion of hydrazine which was considered 


in 
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References [3] and [11], was best represented by n = 2, which is also 
generally true for other monopropellanrs. In this case, a proper 
match of terms can be achieved by setting a' - 0, and balancing 
convection against diffusion. The resulting asymptotic equation 
is : 


d__ 

dr 



ra 

o 


T (1) 


0 


(2-72) 


with the outer boundary condition =0, r = “ from Equation 

(2-69). Therefore, unlike the strand combustion case [3, 11] the 

asymptotic solution for droplets for n = 2 does noc involve the 

reaction rate terra. This occurs since the area for heat conduction 
2 

increases as r rather than remaining constant, increasing the 
significance of heat conduction at large distances from the liquid 


surface. 


The solution of Equation (2-72) satisfying the outer boundary 
condition yields: 


T = 1 - G[1 - exp(-m /r)] (2-73) 

o 0 

at lar^^e r, where the small parameter £ has been absorbed into the 
unknown constant, C. The constant C must be determined by matching 
Equation (2-73) with the numerical solution of Equation (2-56) in 
conjunction with the appropriate boundary conditions at the liquid 
surface. 

Given A, the solution procedure involves guessing values for 

dT 

C and m . Starting values for the numerical solution T , 

o ° I o dr 

are determined from Equation (2-73) at some large but finite value 
of r. Equation (2-56) is then integrated to the droplet surface, r = 1, 
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using the fourth order Hamming predietor-eorrector method. At r = 1 
Equations (2-67) and (2-68) must be satisfied, determining m^ and C 
uniquely. Since Equation (2-56) is nonlinear, a double iteration 
technique must be used. The present procedure employed the Newton- 
Raphson method for this iteration. Once m^ and C are determined, 
the outer boundary condition, Equation ('2-61), is only satisfied 
at the true mathematical infinity, through Equation (2-73), 

for Small A (Small Drop Limit) , For a 
second order reaction, the parameter A in Equation (2-56) is 
proportional to (p^*r^*)^. Therefore, for small droplet sizes, A 
approaches zero and the reaction terms may be neglected in the 
conservation equations. Under these conditions. Equation (2-56) may 
be integrated co yield the following implicit equation for the 
temperature distribution. 

- 1 ) 

i 

The gasification rate m^, may be determined by 

boundary condition, T^ = 1 when r = in Equation (2-74) to yield 

% = i - 'Tos (1 - 1) In [1 - (1 - T^^)/qi (2_ 

The liquid surface temperature, can he determined from Equation 

(2-58). 

Sflution for Large A (Large Drop Limit). At any 
pressure, A increases as the size of the droplet increases. With 
increasing A, the zero-order reaction zone moves closer to the 
droplet, reducing the influence of curvature on the solution. This 
IS also accompanied hy an increase in , although remains 


T + q - 1 

t - 9 -~j 

applying the outer 


unchanged. Reference to the liquid phase solution. Equation (2-63), 
indicates a rapid decay of temperature as one moves into the liquid 
phase, reducing the effect of curvature in this region as well. 

The above discussion indicates that at large A, the process 
approaches a condition where temperature variation and reaction 
effects are confined to a thin region surrounding the liquid surface, 
with little effect of curvature in the spherical coordinate system. 
Under these conditions, the characteristics of droplet combustion 
and strand combustion approach one anotlier and the results of 
References [3] and [11] can be applied directly to the droplet 
combustion problem. Therefore, it is only necessary to convert the 
earlier zero- order results to the presenr notation. 

In References [3] and [11], the expression for the liquid 
surface temperature is equivalent to Equation (2-68) and no change is 
required. Under the present assumptions, particularly unity Lewis 
number, the curvaUnear coordinate system has no influence on the 
surface temperature. This is also true in cases where the Lewis 
number is not unity at the high acti-'?ation energy limit where the 
reaction zone is collapsed to a flame sheet [29], 

In the case of strand combustion, the burning rate eigenvalue 

is given. in terms of a parameter A as follows- 

s 


A 

s 


B* T * A, * 

“O =0 


(p * V *') c * 

000 P 


"o" 1 


M*R*T '* 

OOQ 



(2-76) 


At the large droplet limit, with the combustion zone close to the 
surface, 
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2 

m* = r* 0 *v * 
o s “G «o 


(2~77) 


Therefore, is given by 

*0 (2-78) 

Knowledge of the droplet size end the oheinical kinetic parameters 

provides A; the quantity is given as the eigenvalue of the zero 

order solution in References [3} and [11], completing the 

determination of m „ 
o 

2 . 4 First-Order Problem 

2»4«1 Equations ■ htilizing the general zero-order results, 
the first order equations reduce to the following « 

Liquid Phase r < 1 


= constant 


(2-79) 


dr 


2 *^ ^1 
' dr 

Gas Phase 


m dT, . 2 m,, dT 

Q ■ 1 _ iwr ^ ^ fl p 


dr 6^ 1 dr 

r > 1 


(2-80) 


2 2 2 
T ~r~ - iO)r T = - iwr T 
o dr 1 Q 


d 2 ‘^®1 

^ r T 
d r o dr 


— m “■~S — — = — r 


o dr 




dr 


G dr 

r / 


o 
m 

- m — — + -- — 
o dr dr 


Y - 1 


Y 


101 


r^ T ° 


1 dr 


+ r q w^ 


T r8 - T 

ET 
o 


ri “ 

0 - T 

A G O 


(2-81) 

(2-82) 


dT 


m 


1 dr 


o 2 
-r : 


/Y - 1\ . 2 

^ iw - r q w n 
Y ^ o 


T 

o 


r iu) T, 


(2-83) 
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Boundary Conditions 
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r = 0 = Q 


(2-84) 


r = 1 


m_- = m. T = T, , 

fl 1 Is ls+ 


(2-85) 



- TO. [(l-g)T^ 4- il - TO (l-g)T (2-86) 
i o o i 


= TO C'^, - T.) + TO, (1 - T - q) 
oil 1 o 


0^ - = a q exp(-L^/i^) 


)[^-J 

° It ^ J 


(2-87) 

(2-88) 


Far from the droplet, 0^ and must approach a constant value. 


6l = Ti . 


(2-89) 


The value of depends upon the case under consideration. For an 
unsteady gas phase, since 6^ - large r, the constant can be 

determined by considering the behavior of Equation (2-82) at large 
r. In this case, = (Y-1)/'Y which is the form for isentropic flow. 
In the ease of a quasi-steady gas phase, with transient 
effects still important in the liquid phase, is the amplitude of 
the ambient temperature variation due to the varying energy content 
of the liquid at the surface of the droplet (resulting from transient 
energy storage in the liquid phase). In this case, becomes a 
second eigenvalue in the solution to be determined along with m^. At 
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the limit where the liquid phase also becomes quasi— steady , the 
value of goes to zero, since in the absence of transient energy 
storage in the liquid phase. Equation (2-62) must be satisfied. 

The equations for the first order are linear. Therefore, 
solutions can be combined to eliminate the need for iteration in 
order to determine and K^. Since the limits of integration extend 
to infinity, asymptotic gas phase solutions are employed so that 
numerical computations can be confined to a finite region. Numerical 
difficulties are encountered when integrating the liquid phase 
solution to r = Q, therefore an asymptotic analysis is employed in 
this region as well. 

^'^•2 liquid Phase Solution . Similar to the zero-order case, 
the liquid phase solution is completed by applying the boundary 
condition T^ = T^^ at r = 1, .:ilthough T^^^ is not generally known 
until the gas phase solution is complete. This procedure allows 
the liquid phase solution to be decoupled from the gas phase. 

In order to simplify the subsequent discussion. Equation (2-80) 
may be rewritten as: 


where 






(2-90) 


L = - - . g I- H 

^^2 ^ ^1 dr ^ ®2 


(2-91) 


is a linear operator, and 


H 

1 r 


ra 


PfS^r 




XO) 

6 . 


(2-92) 


(2-93) 
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= (T - T.) 
3 os X 


tfi 

o 


CPjS^r 


2 ^2 


exp 



The boundary conditions on Equation (2-90) are 


(2-94) 


r = 0, 



r - 1, 


T - T. 

1 Is 


(2-95) 


The quantity ih^ is an eigenvalue to be determined from che complete 
gas and liquid phase solutions. Exploiting the linearity of Equation 
(2-90), let 




(2-96) 


Equation (2—96) can be divided into homogeneous and particular parts 





A 



T 


m 


T^* = T^p^’c + (2-97) 

A 

where and are constants to be determined from the boundary 
conditions , The differential equations then become 


= 0 ; = H 3 (2-98) 

A 

The values of and are selected to match the boundary condition 
at ?: = 1 , to yield 


y\ 



= T. /T 
Is' 


IHs- 


K * = - T */T 
T -"iPs- ^^IHs- 


(2-99) 


In order to match the liquxd and gas phase solutions, the 
derivative cf at the liquid surface is required. Combining 
Equations (2—96), (2—97) and (2—99) this quantity is 
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f dr 


s- 


= H, T. + 

4 Is . 5 fl 


(2-100) 


where 


^4 =f 


IHs- 


dl 

IE 

dr 

s- 


«5 " ^fl 


ira.. - laft 


( 2 - 101 ) 


( 2 - 102 ) 


The behavior of Equation (2-98) at small r yields the 
following asymptotic forms for and as r 0 




(2-103) 


T * = : , , 

IP "os 

f f 


(T - T.) 


1 1 1 I 


i-i 

r 


(2-104) 


Equations (2-103) and (2-104) formally satisfy the boundary condition 


^ 0, as r 0, The arbitrary constant in the homogeneous 

yS 

equation has been absorbed into and 

The solution procedure involves employing Equations (2-103) 

and ^2-104) to generate starting values for the numerical integration 

of Equations (2-98) at some small but non-zero value of r. The 

equations are integrated to the liquid surface, allowing the evaluation 

of and from Equations (2-101) and (2-102) , Substitution of 

these quantities into Equation (2-100) provides the needed expression 

in the gas phase boundary conditions . Once T, and ra.., are determined 

Is f 1 

from the gas phase solution, Equations (2-96), (2-97), and (2-98) 
allow complete determination of in the liquid phase. 

2,4,3 Quasi-Steady Gas Phase Solution , A number of 
simplifications are available for the analysis when the gas phase 
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can be assumed to be quasi-steady. This limit arises since 
frequency effects appear as w in the gas phase and w/6^ in the 
liquid phase ^ At atmospheric pressure, for hydrazine, <5^”^ has 
the value 2071, and under these conditions, transient effects in the 
liquid become important at frequencies substantially below the 
frequencies where gas phase transients are encountered. Since 6^ 
is proportional to the mean pressure, the gap between the two 
characteristic frequencies progressively becomes smaller as the 
pressure increases , 


At the quasi-steady gas phase limit, m. = , a constant 

1 fi’ 

For adiabatic combustion, consideration ox the overall energy 
balance and integration of Equation (2-82) shows that 0, 


V a 


constant as well. 

Equation (! 

2-^83) can 

then be simplified 

to yeild 


L = 

'^3 \ ' 

^4 "'l 

-f A^ 

(2-105) 

where 

£» = - 

"2 + \ 

dr 

dr 

"2 

(2-106) 

and 

A = 2 ! 

dT 

^ 4. 

il_ 

ifl 

Q 

(2-107) 


Af Z 

T dr 
. 0 


r T 

0 

1 ‘^^^o 2 

dT w q 

G + 0 

'(5 - n) 


n , E 1 

(2-108) 

2 T . 2 rT 

o dr 0 

dr T 

0 

T 

L 0 

" TT 

- T ) ^ 2 ! 

o- T 1 

0 •' 


q n w 
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q w n 

A, = - -- g 

5 T 


( 2 - 111 ) 


The boundary conditions at r = 1 become 


dr ^2 "^Is " ^3 \ 


(2-112) 


dT. 


' ' "L ■" ^2 ’^1 ■ ^3 “l + ”l 


1 1 dr 


(2-113) 


■'l ■ >^1 ^ «2 


(2-114) 


where 


^1 = 


T 


os 


(2-115) 


dT \ 

dr - «4 

(2-116) 

= (1-$)T + L + 

os 5 

(2-117) 

F_ ^ - T 
1 os 

(2-118) 

/dT 1 


(2-119) 

' »s- 

= m 
3 o 

(2-120) 

F. = (1 - T - q) 
4 os 

(2-121) 

a q L- 

+ exp(-L /T ) 

F£ ^ X O S 

(2-122) 


o 


= a q exp(-L^/T^^) 


The outer boundary condition is 


r T^ ^ 


Again exploiting linearity, define 


(2-123) 
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Ti = + th^ (2-124) 

Separating the functions in Equation (2-124) into homogeneous and 
particular parts 






✓v 


+ K T 


IH 


T * = 
1 


T * + K* T 
IP ^IH 


where 


Ti . lu, + K T^j, 

(2-125) 

i'{V - 0 

i'iii,) - A^ 


£'{Ti/} - 

i’CTip} = 

(2-126) 


The parameters K, K* and K are constants which can be selected to 


satisfy any of the boundary conditions. Since the outer boundary 
condition can be handled by the asymptotic solution of Equations 

s Equation (2-114) was chosen to determine the values of the 
K's. This yields 


K = 

\G^ IPs+j 

/T 

' ' lHs+ 

(2-127) 


= _T * /T 

IPs '^^1HS+ 

(2-128) 

K = 

^1 - ' 
[Gi lPs+1 

|^'‘^lHs+ 

(2-129) 


The remaining two boundary conditions at r = 1, Equations (2-112) 
and (2-^113), can be employed to determine the values of and K^, 


The result is 
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\ ^'^11 '^23 ■ '^21 '^13^'^ ^'^11 “^22 ■ '^12 ‘^21^ (2-130) 


™1 ~ ^'^13 '^22 ~ '^12 “^23^ '^^■^11 *^22 " '^12 ‘^21^ (2-131) 


where 


dT *\ 

j = D 1 — ^ + D T * - D 

11 1 I dr 2 ls+ 3 

St 


dT 

h2 - h ldT-|,H. “2 h.+* 


h3 ■ - ®1 - »2 h. 


21 


1 1 L 


+ ^2 ^is" - ^4 


dTi*\ 

‘^22 ^1 I dr J ^2 ''■’ls+ ^3 

' s+ 


dT 

hs = h \l^l^ - ^2 h.+ 


(2-132) 


The derivatives in these expressions are obtained by differentiating 
Equations (2-125) noting that the K’s are eonstantSc 

In order to complete the solution, the asymptotic forms of 
Equations (2-126) at large r must be solved. Proceeding in the 
same manner as the zero-order case, the result is 


r ->■ CO 


^IH ” ” ^IP* ~ ^ exp(-m^/r) 


Tip = exp(-ih^/r) 


(2-133) 


Equations (2-133) automatically satisfy the outer boundary condition 
of the gas phase solution. 
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The solution is obtained by calculating D^, D^, F^, , 

^31 2 nd from the zero-order solution and the solution of 

the first order liquid phase. The parameters - A^ are also 
available from the zero order solution. Equations (2-133) are used 
Co generate starting values for the differential equations given by 
Equations (2-126) at some large but finite value of r. These 
equations are then integrated to the surface of the droplet. The 


parameters K, K*, and K can then be determined from Equations (2-127) 

/N 

to (2-129), When these parameters are ascertained, T. . T *. 

’ is’ Is * 

T^^, and the derivatives of these quantities can be determined from 
Equations (2-125) and their derivatives. Determining the J.. terms 
from Equations (2-132) then allows and ra^ to be found from Equations 
(2-130) and (2-131), If desired, T^ can then be evaluated as a 
function of r from Equation (2-124) for the gas phase (yielding T ) 
and from Equation (2-96) for the liquid phase. 


2.5 Wet Sulb Limit Equations 

In the early phases of the lifetime of a droplet, and for 
large droplets having strong decomposition effects, substantial 
temperature gradients exist within a drople c [30] , The present 
investigation examines this limit by considering liquid temperatures 
at Che center of the droplet that are lower than the equilibrium 
surface temperature. 

In the later stages of droplet gasification, particularly 
for small droplets with relatively weak monopropellant decomposition 
effects , temperatures throughout the droplet approach the liquid 
surface temperature [30]. This can also occur due to circulation 



1 




/ 
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in the liquid phase; however, this effect is difficult to determine 
quantitatively. When temperatures throughout the droplet are equal 
to the surface temperature, droplet conditions are often referred 
to as combustion at the wet bulb state. In order to cover the complete 
range of possibilities, the present calculations also consider 
response at the wet bulb state. In both cases; however, the overall 
combustion process is assumed to be adiabatic. 

The distinction in the determination of property values for 
the temperature gradient and wet bulb cases may be seen as follows, 

For adiabatic combustion, the flame temperature is given by the 
dimet-sional form of Equation (2-62) 


W' " ^ ^f* V ■ (2-134) 

With temperature gradients present, can be specified (as the 

liquid injection temperature, for example) and the flame temperature 

can then be obtained from Equation (2-134), The other dimensionless 

variables are determined from Equations (2-22) , and the liquid 

surface temperature from Equation (2-68). 

For wet bulb combustion, liquid temperatures are uniform and 

T^* = 3ud both Tj_* and are unknown. In this case, 

Equation (2-134) and the dimensional form of Equation (2-68) are 

solved simultaneously to yield T.* and T *, followed by the 

1 «o 

determination of the dimensionless variables through Equations 
( 2 - 22 ) . 



CHAPTER III 


STEADY STATE DROPLET RESULTS 


3,1 Introduction 

The physical property values and the chemical kinetic 
parameters used in the calculations are representative of hydrazine 
decomposition. The specific values, summarized in Table 1, are 
identical to those used by Allison and Eaeth [11], 

A second order reaction correlates the strand burning rate 
of hydrazine at pressures above 1 atmosphere [11], For steady state 
conditions, the theoretical burning rate is relatively insensitive 
to variations in the activation energy. However, for oscillatory 
combustion, the activation energy influences both the amplitude 
and phase of the combustion response of strands. Allison and Faeth 
[11] found than the value of E* given in Table 1 gave good 
agreement between rheoreticai and experimental response determinations 
over their test range. This value also agrees with earlier values 
reported for hydrazine decomposition [31], With the activation 
energy given, the value of B* is found by matching the steady strand 
burning rates - 

3,2 Temperature Distribution 

In order to check the accuracy of the zero-order model, 
predicted and measured liquid surface temperatures were compared as 
a function of pressure. The theoretical prediction of liquid surface 
temperatures is independent of geometry under the unity Lewis number 
assumption. Therefore, the comparison is identical tc that made in 
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Table 1 


Properties Used in the Theoretical Model 


Property 

Value 

a*(N/m^) 

3.0078 X 10^^ 

B* (m^/Kg-s) 

8.385 X 10^° 

C^*(KJ/Kg-K) 

3.0961 

Cp* (KJ/Kg-K) 

3.0961 

E*(KJ/Kg-mol) 

111700 

L*(KJ/Kg) 

1715 

L^*(KJ/KG-raol) 

40794 

Le 

1 

M*(Kg/Kg-mol) 

24 

n 

2 

q*(KJ/Kg) 

4-.37 

R* (Kj/Kg-mo].-K) 

8.3144 

T^* (K) 

298 

T *Co) 

1345 

Y 

1.126 

6 

0 

A^*(J/m-^s-‘K) 

0.390 

A„^A(JAn-s^k) 

Pj*(Kg/m^) 

0.176 

1000 
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Reference [11] between measured and predicted surface temperatures 
in the strand combustion case. These results have been replotted 
and are shown in Figure 2. The agreement is good throughour the 
pressure range of the measurements (.51-19.8 atm.). 

The test conditions shown in Figure 2 involved a liquid 
temperature of 298 “K far from the liquid surface. This is 
representative of a non-wet bulb condition, with the liquid temperature 
increasing as the surface is approached. A theoretical calculation 
for the wet bulb state is also presented in Figure 2, where the 
liquid temperature is constant. For the wet bulb case, no energy 
is required to raise the liquid to the vaporization temperature, 
resulting in a higher surface temperature than the non-wet bulb case 
(for a given pressure). 

Steady state calculations were completed over a range of A 
(essentially representing droplet size variations) and pressure 
for both non-wet bulb (with a centerline droplet temperature of 298°K) 
and wet bulb conditions. 

Figures 3 and 4 illustrate liquid and gas phase temperature 
profiles for non-wet bulb conditions at a pressure of lO^N/m^ 
(appro.icimately 10 atm,). Curves are plotted for various values of A, 
When A = 0, there is negligible reaction in the gas phase surrounding 
the droplet, and the process corresponds to simple droplet evaporation, 
in this ease, temperature variations are gradual and the effect of 
cur\7-ature in the spherical coordinate system is important. As A 
(droplet size) increases at a given pressure, the rate of 
gasification also increases. The greater liquid flow rate toward the 
surface of the droplet tends to sweep back the temperature profile 



P/xlO^®(N/m^) 


Figure 2 Theoretical a.nd Experimental Liquid Surface Temperatures as a Function of Pressure 
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in the liquid phase (Figure 3) . The increased reaction rate also 
causes the combustion zone in the gas phase to approach the liquid 
surface as shown in Figure 4. Therefore, as A increases, temperature 
gradients become confined to the region near the liquid surface and 
the effect of curvature is reduced. These conditions correspond to 
the large drop limit, where strand combustion is similar to adiabatic 
droplet combustion. 

Figure 5 shows gas phase temperature distributions for the wet 
bulb case, with remaining conditions similar to the calculations shown 
in Figure 4, For wet-bulb conditions the liquid temperature is 
constant and equal to the surface temperature. These gas phase 
temperature profiles are similar to the non-wet bulb case, becoming 
steeper as the value of A increases , The maj or difference is the 
slightly higher temperature at the liquid surface for wet bulb 
conditions , 

5 2 7 2 

Calculations completed for pressures of 10 N/m and 10 li/m 
(approximately 1 and 100 atm,) gave results that are generally 
similar to those shown in Figures 3-5, 


3.3 Mass Burning Fates 

The mass burning rate of a liquid droplet depends on the size 
of the droplet (at a given pressure) and approaches two asymptotic 
limits: the large drop and small drop limit burning rates. The 

small drop limit represents a droplet vaporizing without reaction, 
and the large drop limit approximates a one-dimensional liquid strand 
with the flame very close to the liquid surface. The intermediate 
drop sizes or transition range constitute the majority of drops that 
are distributed between the two asymptotic limits. 




?5.gure 5 Steady State Gas Phase Temperatures as a Punction of A ah Approximately 
10 Atm, Pressure (Wet Bulh Conditions) 
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Figure 6 shows the steady state mass burning rates as a 

function of A at 1 atmosphere pressure for non-wet bulb conditions. 

The complete steady state solution clearly approaches the large and 

small drop limits. In order to plot m^ as a function of A at the large 

drop limit, values of A^ were obtained from the strand combustion 

calculations [11] , Equation (2-^78) then yields the large drop limit 

mass burning rate, m , as a function of A. The small drop limit 

LBL 

mass burning rate, m , is calculated directly from Equations 

SDL 

(2-68) and (2-75) and is independent of A. 

An approximation to the complete steady state solution is 
shown in Figure 6 . The approximation was obtained by simply adding 
the small and large drop mass burning rates at each value of A, as 
follows ; 


o 

m 

o 


m + m 
°SDL °LDL 


(3-1) 


The approximation of Equation (3-1) yields a maximum error of 
approximately 12% over the entire range of A. Equation (3-1) will 
be used later to evaluate spray combustion characteristics since it 
can be readily integrated to determine droplet lifetimes. 

Data on hydrazine combustion was available from experimental 
work done by Allison and Faeth [9] , This study involved droplet 
burning in a combustion gas under decomposition conditions, for 
varrous droplet sizes, at atmospheric pressure. The data was 
limited to the ambient temperature range 1660-2530“K; therefore, it 
was necessary to extrapolate the measurements to the 1345 ’’K ambient 
temperature of the present study. The values used in the experimencal 
comparxson shown in Figure 6 are summarized in Table 2. The agreement 


P* =10® N/m* 

Tj* = 298K 
T«o »I345K 

O EXPERIMENTAL, REF 9 


COMPLETE SOLUTION y 

approximation— V \ 



'LARGE DROP 
LIMIT 


■SMALL DROP LIMIT 


igure 0 Theorecical and Experimental Steady State Droplet Ga«ifica-i or, R ■ 
« 1 Atn.. P^ssur. (Nop-, tec Bulb lopdnSor 




between the predicted and measured burning rates is seen to be 
reasonably good, lending confidence to the theoretical model. 


Table 2 


Steady Droplet Burning Rate Data (Non-Wet Bulb)^ 


d^* (diameter) 
(mm) 

Awiiio* X 10^ 
(Kg/s) 

1.3 

,70 

6,3 

1,4 

19.1 

135 


Extrapolated to = 1345K from data of Reference [9] ar 

atmospheric pressure. Non-wet bulb conditions, T,* = 298K, 

Burning rate results at atmospheric pressure for wet bulb 
conditions are illustrated in Figure 7, The general behavior is 
similar to the non-wet bulb conditions; the major difference involves 

a slight increase in the burning rate throughout the entire range of 

A. 

Theoretical mass burning rates at pressures of 10 and 100 
atmpsheres were completed and gave results very similar to Figures 6 
and 7. In all cases the transition region between the large and small 
droplet limits falls approximately in the range 10^ < A < 10^. The 
approximation of Equation (3-1) , was found to provide a good 
representarioB of the data (with maximum errors on the order of 10%) 
over tue entire range of the calculations. 
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In order to employ the approximate burning rate expression, 

it is necessary to have values for m and A , These quantities, 

°SDL ® 

along with surface temperature predictions , are summarized in Table 
3, over a range of pressures, for both wet bulb and non-wet bulb 
conditions , 


Table 3 


Summary of Steady State Calculations 



p 

o 

T. 

1 

T 

os 

T *(K) 
os 

T *(K) 

ih 

°SDL 

A 

s 


Non 

Wet Bulb 

Conditions 





10^ 

.221 

.278 

374 

1345 

.544 

1.432 X 10^ 

10^ 

.221 

.335 

450 

1345 

.509 

1.432 X 10^ 

10^ 

.221 

.418 

562 

1345 

.434 

1,432 X 10^ 

Wet 

Bulb Conditions 





10^ 

.266 

,266 

376 

1417 

,603 

6,205 X 10^ 

10^ 

,306 

,306 

457 

1498 

.629 

2,660 X 10^ 

10^ 

.359 

.359 

582 

1623 

.668 

8,162 X 10^ 


With : 

Increasing 

pressure , 

the surface 

temperature 

increases in 


all cases. As noted earlier, the wet bulb surface temperature 
slightly exceeds the non-wet bulb surface temperature at each pressure. 
The nondimens ional surface temperatures are lower for the wet bulb 
condition at a given pressure even though the dimensional surface 
temperatures are higher. This is due to the fact that at the wet 
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bulb condition, the flame temperat' 4 i*es and ambient temperatures 

are also higher. At the non-wet bulb limit, for fixed 
centerline temperature, the ambient temperature is independent of 
pressure , 

I'/hile A^ is relatively constant with pressure at non-wet 
bulb conditions, A^ decreases with increasing pressure for the wet 
bulb case, significantly increasing the mass burning rates for the 
larger drops at a given value of A. This effect is due to increased 
reaction rates for the higher flame temperatures of wet bulb combustion, 
at elevated pressures. 

At non-wet bulb conditions , m decreases with increasing 

°SDL 

pressure due to the reduction in the temperature difference between 
the ambient gas and the liquid surface, For the wet bulb case, the 
increased ambient temperature compensates for this effect and thera 

is a slight increase in m with increasing pressure. 

°SDL 

Using the properties of Table 1, the droplet siae range 

4 9 

corresponding to the transition region (taken to be 10 < A < 10 ) 

is shown in Table 4, As the pressure increases, the droplet size 
for the onset of the large drop limit is reduced, Notably, the size 
range of technological importance for actual combustors falls 
largely in the transition region. At higher pressures, however, a 
greater percentage of the droplets present in a spray can be 
represented by large drop limit results. 



Droplet Size 


P * X 10^ 
o 

(N/m^) 

1 

10 

100 

3 . 4*9 

Corresponding to 10 < A < 10 
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CHAPTER IV 


EIRST-ORDER RESULTS 


4«1 Introduction 

At very low frequencies , the liquid and gas phases respond 
rapidly to changes in the forcing function and a total quasi-steady 
analysis is sufficient c As the oscillation frequency is increased, 
the liquid response is not rapid enough as the characteristic 
frequency of the liquid phase thermal wave is approached, requiring 
an unsteady analysis of the liquid phases As the frequency approaches 
the characteristic frequency of the gas phase thermal wave, both the 
liquid, and gas phases are unsteady. At extremely high frequencies, 
both the liquid and gas phases are unable to respond at all, and a 
steady state approximation can be assumed- 

In order to analyze response, the burning rate response 
function is utilized [1]« In the present notation this quantity is 

= Re{m^/ifi^} (4-1) 

where Re denotes the real part; i.e,, that portion of the burning 
rate oscillation that is in-phase with the pressure oscillation, Eor 
instability, P_ must be positive and of order unity (the exact value 
depends upon the degree of damping present) at a point in the 
combustion chamber where the pressure is varying. It must be 
understood that a droplet with ^ 1 will not in itself cause 
instability, but rather the sum of the responses of all the droplets 
within the combustion chamber must be greater than the a'''’ailable 
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damping, in order for the combustion oscillation to 


grow. 


4 . 2 Small Drop, Low Freq.uency Limit 

This section considers behavior at low frequencies for the 
small droplet or evaporative limit. Figure 8 illustrates both the 
wet and non-wet bulb response function, and the magnitude of 
the surface temperature oscillation, plotted as a function of 

pressure. The frequency for this plot is zero so that both the 
gas and liquid phases are quasi-steady. (The present analysis is not 
formally valid at w = o, since the mean droplet size varies during an 
oscillation at this condition; therefore, the plot should be taken 
as representative of conditions where it is still valid to assume 
that both gas and: liquid phases are quasi-steady, but with sufficiently 
hrgh frequency so that changes in mean size can he ignored during a 
period of oscillation.) 

■A.t this condicion, the response function is negative, with 
a relatively small magnitude, and with increasing pressure, the 
response function becomes more negative. This is in contrast to 
results for biprcpellant comhustion where the response function 
generally approaches zero or a positive limit for small u. 

At this low frequency limit, the ambient temperature is 
constant since there is no transient energy storage in either phase, 
llie surface tt jperarure varies in phase with the pressure oscillation, 
and its magnitude increases as the pressure increases. This follows 
primarily from the form of the vapor pressure curve. At the quasi- 
steady limit 
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T 


Is 


ps 



(4-2) 


therefore, since increases with increasing mean pressure, 
increases accordingly . 


At the same limit, using the results for T. 


IIL - 


os 


1 ' L /.1 + 


(4-3) 


This approximate formula indicates negative response at the quasi- 
steady limit. With increasing pressure, increases, providing a 
larger negative response. 


Physically, the negative response at the totally quasi— steady 
limit is caused by the fact that the gasification rate decreases as 
the surface temperature increases (for a constant ambient temperature) . 
Therefore, pressures above the mean result in reduced gasification 
rates, and pressures below the mean result in increased gasification 
rates -yielding a negative response. At higher frequencies, energy 
can be temporarily stored in the liquid or gas phases, causing a 
fluctuation In the ambient temperature which can compensate for the 
Surface temperature fluctuation. This modifies the response for 
small droplets as will be discussed later. 


4.3 Large Drop Limit 

Moving from the evaporative limit (small drop limic) to the 
large drop limit, much greater response is observed for monpropellant 
droplets. Figures 9—11 illustra.te the non— wet bulb response as a 
function of frequency (dimensionless frequency for strand combustion 
from Reference [11]) for the large droplet limit at mean pressures of 






l?±gure 10 Response as a Runction of Frequency at the Large Droplet Limit at 
Approximately 10 Atm„ Pressure (Non-Wet Bulb Conditions) 
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1, 10 and 100 atmospheres. Figures 12-14 show similar results at 
wet bulh conditions. In addition to the quasi-steady gas phase 
approximation, results for the completely unsteady gas phase are also 
shown. The non-wet hulh results are replotted from the calculations 
of Reference [11], v^hile the wet hulb results were calculated by the 
present author. 

The conversion of r^ and r^^, the strand burning rate and its 


perturbation, is straight^forv<rard at the large drop limit, namely 



Re {r^} 


(4-4) 


since r^ = 1. The dimensionless frequency for the strand combustion 


case, is related to the present dimensionless frequency by 


Wg = (A^/A)w (4-5) 

In all cases, the response approaches unity at low frequencies, 
where both liquid and gas phases are quasi-steady. This follows from 
the fact that the burning rate is proportional to pressure, under 
steady conditions, for a second-order reaction [11]. 

With increasing frequency, at non-wet bulb conditions, a 
peak is observed in the response plot at frequencies near the 
characteristic frequency of the liquid phase. Beyond this peak, the 
quasi-steady gas phase solution gives a continuously declining 
response, The peak is absent in the wet bulb case, with the response 
showing a noticable decline at the liquid phase characteristic 
frequency. This increasing contrast between the two cases will be 


discussed later in more detail. 


* ® 



Figure 12 Response as a Function of Frequency at the Large Droplet Limit at 
Approximately 1 Atm=, Pressure (Wet Bulb Conditions) 




i 






I'igure 14 Response as a Junction of Frequency at ttie Large Droplet Limit at 
A.ppr.::-<iinately 100 Atm. Pressure (Wet Bulb Conditions) 
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At higher frequencies, the analysis allowing for transient 
gas phase effects begins to diverge from the quasi-steady gas phase 
solution. Substantial differences between the two models are 
encountered for values of on the order of lO”^, which represents 
frequencies near the characteristic frequency of the gas phase. At 
dimensionless frequencies on the order of unity, a second response 
peak is observed, in all cases, for the complete unsteady gas phase 
solution. This peak is clearly associated with gas phase transient 
effects and yields maximum values of on the order of unity. 

With increasing pressure, the frequencies of the two response peaks 
approach one another due io the increase in gas density. 

Table 5 lists the frequencies from Figures 9-14 where large 
droplets have a positive response, the liquid thermal wave 
characteristic frequencies, and the gas thermal wave characteristic 
frequencies. In considering this table, it should be recalled that 
acoustic instability in rocket engine combustion chambers is 
generally associated with the frequency range 500-30,000 Hz [1] 

At pressures below 10 atmospheres, response peaks due tc gas phase 
transient effects largely fall within tills range. At pressures above 
10 atmospheres, response peaks due to gas phase transient effects 
decline in importance since they are associated with very high 
frequencies- In this pressure range, liquid phase phenomena are 
more significant, providing response peaks in the critical frequency- 
range, Since the pressure range usually encountered in rocket 
engines is above 10 atmospheres, gas phase transient effects do not 
appear to be a major factor in causing instability for hydrazine 
fueled engines. Based on this finding, further analysis in the 
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Table 5 
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Critical Frequency Ranges for Hydrazine Strand Combustion 


(Non-Wet Bulb Conditions ) 


Mean 

Pressure 

(N/m^) 

Liquid Thermal 
Wave 

Characteristic 

Frequency 

(Hz) 

Liquid 

Transient 

Range 

(Hz) 

Gas Thermal 
Wave 

Char acteris tic 
Frequency 
(Hz) 

Gas 

Transient 

Range 

(Hz) 

10’ 

0.02 

< 1.8 

210 

270-1,900 

10*^ 

2.1 

< 84^ 

2,100 

2,100-170,000 

10^ 

210 

< 1,800^ 

21,000 

27,000 


(Wet Bulb 

Conaitions) 




Mean 

Pressure 

(N/m^) 

Liquid Thermal 
Wave 

Characteristic 

Frequency 

(Hz) 

Liquid 

Transient 

Range 

(Hz) 

Gas Thermal 
Wave 

Characteris tic 
Frequency 
(Hz) 

Gas 

Transient 

Range 

(Hz) 

10^ 

0.02 

< 2 

210 

330-1500 

10^ 

1,8 

< 130 

2,100 

2,100-18,000 

10^ 

180 

< 6700'^ 

21,000 

6,700^-190.000 


liquid transient range ends where transient gas curve intersects 

R = 0. 

r 

b„ . 

Gas transient range begins where second gas transient peak intersects 

P = o. 
r 

liquid transient range ends where gas transient range begins ^ at 
low point in gas transient curve at approximately P = 0.14. 
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transition region, between the large and small droplet regimes, was 
limited to the case oj quasi-steady gas phase effects. 

The magnitude of the unsteady gas phase response peaks 
approach unity at elevated pressures; however, these response peaks 
are smaller In magnitude than the transient liquid phase response 
peak for the non-wet bulb ease, Also, the unsteady gas phase peak 
otcurs near or beyond the maximum probable frequencies that cause 
instability. Table 6 lists the frequencies where the unsteady liquid 
and gas phase response peaks occur. At the higher pressures, the 

unsteady gas phase response peak is outside the critical Instability 
frequency range. 

Table 6 


Frequencies of the Unsteady Gas and Liquid Phase Response 


Peaks 



(Non-Wer Bulb Condition) 
10^ 

10 *^ 

10 ^ 


(Wet Bulb Condition) 
10^ 

10 *^ 

10 ^ 


0.U4 

17 

1100 ^ 


Corresponds with first gas phase transient analysis peak. 


528 

11,000 

110,000 


670 

6,700 

67,000 
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4.4 Complete Solution 

Figures 15-17 show the quasi-steady gas phase non-wet bulb 
response as a function cf frequency at various drop sizes, ranging 
from the small drop to the large drop limit. As A increases, the 
response increases, caused by the combustion zone moving closer to 
the liquid surface. At an A of approximately 10^, a peak in the 
response curve develops, increasing in magnitude and frequency as 
the drop size increases, until the large drop approximates the one- 
dimensional strand near the value of A = 10^^, As the pressure 
increases, the peaks increase in magnitude and move to higher 
frequencies, it should be noted that when A = 10 , transient gas 
phase effects begin to appear at approximately (jJ = 6 98 x 10^ „ This 
corresponds to - 10~^, 

The quasi-steady wet bulb responses are plotted in Figures 
18—20, as a function of frequency for the complete range of drop 
sizes. As in the non-wet bulb case, the response increases with 
increasing A; however, no response peaks are observed. Droplet 
combustion response is well approximated by the one-dimensional strand 
at values of A greater than lO^'^. 

In order to better understand the difference in the response for 
the wet and non— wet bulb cases, an investigation of the first-order 
"liquid phase temperature distribution was undertaken. 

Figure 21 shows the non-wet bulb perturbation surface 
temperature, magnitude and phase angle, as a Eunction of frequency 
f ot a pressure of 10 atmospheres and A == 10 ^ , At low frequencies 
the surface temperature, is in-phase with the pressure 

oscillation, l^hen the frequency increases, a peak occurs in the 
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Figure 20 Quasi-i:teady Gas Phase Response as a Function of Frequency and A 
at Approximately 100 Atm. Pressure (Wet Bulb Conditions) 
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Figure 21 Perturbation. Surface Temperature Magnitude and Phase Angle as a 
Function of Frequency and Pressure (Non-Wet Bulb Conditions) 
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magnitude of T^^, whieh corresponds with the response peak. It is 

also observed, that when the peak does occur, T leads the forcing 

xs 

function until the peak is reached, while at higher frequencies, T 
lags the pressure oscillation. 

Figure 22 is a plot of the first-order liquid phase temperature 

distribution near the drop surface for a sinusoidal oscillation at a 

nondimensional frequency of 1,6, a pressure of ID atmospheres, and 

A = IQ^^. This plot further illustrates the phase lead. A 

frequency of 1,6 corresponds to the maximum phase lead in Figure 21 

for the perturbation surface temperature at 10 atmospheres. As the 

pressure wave oscillates at o = 1,6, T, also oscillates at 

Is 

approximately 4 degrees ahead of the pressure , is unable to 

stay in-phase with the pressure oscillation at this frequency because 
the large values of induce steep temperature gradients near the 
liquid surface which act as heat sinks . The increased heat loss 
from the surface due to this gradient causes the temperature peak to 
precede the pressure peak for this range of frequencies (and 
conversely at che minimum pressure and temperature regions). 

Figure 23 shows the magnitude and phase angle of the wet 
hulb perturbation surface temperature as a function of frequency at 
10 atmospheres for A = 10^*^, Ic Is observed that T^^ never leads the 
pressure oscillation in this case, in contrast to non-wet bulb 
conditions (Figure 21) . 

In order to better understand this phenomena, a plot similar 
to Figure 22 is introduced. Figure 24 is an illustration of the 
liquid: phase first-order temperature distribution near the drop 
surface at w = 50 , a pressure of 10 atmospheres, and A = 10^*^. In 
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Figure 22 Perturbation Liquid Phase Temperature Distribution for a Sinusoidal 
Pressure Oscillation and a Tis Phase Lead of Approximately 4 Degrees 
(Non-Wet Bulb Conditions) 
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this case lags the pressure oscillation by approximately 5.5 
degrees. The wet bulb case differs from the non-wet bulb condition 
because of the constant zero-order mean liquid phase temperature, 
wnich reduces the steepness of the first-order temperature gradients 
in the liquid phase near the surface. 

The only two types of monopropellant droplets which have been 
considered, are the wet bulb droplet with a zero mean temperature 
gradient (constant mean temperature) and the non-wet bulb droplet with 
a steep mean temperature gradient. For the parameters listed in 
Table 1, a response peak is observed in the non-wet bulb case, while 
the opposite is true for the wet bulb case. Under realistic 
conditions, the droplets will have temperature gradients ranging 
between the wet and non-wet bulb limits. 

From the present study it can be concluded chat a response 
peak will only occur when the mean temperature gradient near the 
droplet surface has sufficient steepness so that enough energy will 
flow from the surface towards the center of the drop, causing the 
perturbation surface temperature to lead the pressure. When the 
mean temperature gradient is not steep enough, the perturbation 
surface temperature will never lead the pressure, A-lso, as the 
droplet size decreases while it is reacting, the liquid temperature 

gradient will tend to approach the wet bulb limit; thus, eliminating 
the response peak. 

Once the response peak begins to appear, leads the pressure 

wave. As the frequency is increased, the non-wet bulb temperature 
gradient loses its influence on the droplet surface, enabling the 
\s to begin to decrease, l^lren the frequency is further 
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increased, the first-order surface temperature and the pressure wave 
are again in-phase; thus, the burning rate is maximized and the 
response peak is reached. Any additional increase in frequency 
will decrease the response as the lag in surface temperature 
continues to increase. 


CHAPTER V 


SPRAY COMBUSTION MODEL 


As discussed earlier, a eomlDustion chamber will become unstable 
when the total response of all the droplets in the spray exceeds the 
available damping , Using the droplet response results of the 
urevxous chapter, an analysis of spray response will not be undertaken. 


5-1 Description of the General Model 

For the spray combustion model, the spray is assumed to be 
relatively eonfined, with the pressure, field identical for all parts 
of the spray. 


The open-loop response function for an individual drop is 
defined on page 130 of Reference [1] as follows: 


P = R 


J P* 


(5-1) 


P*' m* 

where P** and m*' are the instantaneous perturbations of pressure and 
burning rate, and P* and in* are the respective steady state values, 
For a monodisperse spray, Equation (5-1) can be extended to 
compute the total response of a spray consisting of numerous drops 
of varying sizes as follows; 


A-.*' w.. 

P = R / ^ 

rt eVP*’ 




(5-2) 


where ih^* and m^*’ are the respective total steady state and 
perturbation mass burning rates of all the drops in the spray, 
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The total number of droplets in a spray, N, is ajiiul tc the 

droplet injection rate multiplied by the droplet lifetLme, The 

number of droplets in a size increment dr * is Nfdr *. where f is a 

s s'"' 

function of r^*, pressure, and the combustion charaGteristics . The 
fraction of droplets in a given interval dr * is fdr *. 

3 s’ 

With these definitions, the total steady state mass burning 
rate for all drops in the spray is 


DO 

V ■ V* = / 


Nf4Hm * dr * 
o s 


(5-3) 


where ih^* is the steady state mass burnihg rate per unit solid 
angle and is a function of r ^ at a given pressure) 

The total mass burning rate perturbation is 


00 

m * = / Nf4nth * dr * 
1 s 


(5-4) 


where 


. . P* m* ' 

in ^ = — ■ '• ■ — 

- p*, 


m * = 5fj *' 

1^ P*' T 


(5-5) 


and m^* is the perturbation mass burning rate per unit solid angle 
and is a function of r^*, frequency, and pressure for the quasi- 
steady gas phase ease. 

Substituting Equations (5-3) and (5-4) into Equation (5-2) 
yields the total spray response, P^^, for a monodisperse spray 

.'I 


Nf 4Eih, A dr * 
1 s 


P , = R 
rt e 1 “ 




Nf4IIm * dr * 
o s 


(5-6) 
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Introducing the dimensionless radius 


R = r */r * 

s s so 


(5-7) 


where is the radius of the droplet initially injected into the 

combustion chamber; recalling the definition of the dimensionless 
mass burning rate from Equation (2-22); and noting 


Re {m. } = P m 
1 r 0 


(5-8) 


Equation (5-6) becomes 




(5-9) 


Equation (5-9) can be solved numerically by dividing the spray 

distribution into finite increments of AR and summing the weighted 

s 


responses of all the size groups . 


When a polydisperse spray (more than one initial size droplet) 
is injected into the combustion chamber, the overall spray can be 
thought of as many separate spray sets being injected into the 


chamber, where each spray set is a monodisperse spray. 


The number of droplets from each spray set in an increment 
dtg* is Nfdr^*, where N is the total number of droplets in each spray 
set and is related to the number of drops injected per unit time per 
set, and f has the same definition as for a monodisperse spray. 

The number of drops injected per unit time in an initial size 
interval dr^^ * is Igdr^^*, where I is the total number of drops 


1 
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injected per unit time, and g is a function of r *, pressure, and 

s o 

the injector. The fraction of droplets injected per unit time in 

a given interval dr * is gdr *. 

° so ° so 

Two new dimensionless radii are defined as follows: 


s 


s r * 

ref so 

ref 


* 


R 


so 


so 


ref 


so 


ref 


0 < R, < 1 

’ref 


0 < R 
— so 


C 1 


ref 


(5-10) 


where r * is the largest drop injected into the combustion 
®®ref 
chamber. 

Substituting Equation (5-10) into Equation (5-9) and taking 
into account the polydisperse spray, gives the following result 


P 

rt 


1 

f ° 

- 1 

f NfP ra R dR 

/ r o s _ s ^ 

J ref ref 

t-'q J 

dR 

ref 

1 

r 1 


f g 

I Nfih R dR 

dR 

J 

J ° ®ref ®ref 

so „ 
ref 

o 

1 

- o 





(5-11) 


The polydisperse spray can be solved numerically by separating 
the spray into drop size groups and individual spray sets, by a 
simple extension of the method used for monodisperse sprays. 

The remaining assumptions are similar to those employed in the 
response analysis. The effect of droplet heat-up is ignored, with 
temperature gradients assumed to persist in the liquid phase. The 
effect of convection is also neglected so that specific combustion 
chamber configurations need not be considered at this point. 


BfiPBCTODCmLOT 

OBTGINAL PA-GE IS 

The distribution of droplet sizes in the combustion chamber 
depends upon the range of droplet sizes injected and the steady 
Gombustion characteristics of each size, figure 25 outlines the 
process, schematically. Each spray set has a total lifetime, t^^*. 
The number of droplets m the chamber originating from a given 
initial size is equal to the number of droplets of that size 

injected per unit time, times the lifetime of that size. The number 
of drops per size group, resulting from a given initial size, is 

given by the fraction of the lifetime spent within any subsidiary 
size group. 

When more than one size droplet is injected, the additional 
spray sets can contribute to droplet size groups of the primary spray 
set as indicated in Figure 25. In this case, the size distribution 

of the droplets within the spray must be obtained by summing over all 
drop size groups , 

^ hife History Analysis 

As indicated earlier, the time variation of droplet size must 
be known in order to determine the spectrum of droplet sizes in a 
spray from the injected size distribution. This can be determined 
by considering the mass balance of the droplet 

* 


/,Tr +2 dr * 

s = - 4M* 
® f dP^ • 




where 


A/Cp* 


(5-13) 


The dimensionless steady state mass flow rate per unit solid 
angle for any given sized droplet, can be approximated (as explained 
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in Chapter III ) by the sum of the small drop limit and the large drop 
limit mass flow rates as follows; 


m = m 
o o. 


+ ih 


(5-lA) 


'SDL ^LDL 

The large drop limit mass flow rate is a function of from 

Equation (2-78) and can be transformed into 


m 


LDL 



[V 1 

*^s ' 

r *1 
so / 


(5-15) 


where r^^* is the initial drop radius, A is the value of A for r * 


So 


one- 


and the combustion chamber pressure, A is obtained from the 

s 

dimensional strand analysis, and A is calculated from Equation 

SDL 

(2-75). Thus, the mass flow rate of a given sized droplet can be 
approximated by 


m = ra + (R ) 

® °SDL ' 1 ® 


(5-16) 


In order to put Equation (5-^12) into a generalized form, R and 

s 

(5-17) 


T = 


t*A * 
0 ^ 


r * p^*G * 
so f P 

are introduced. 

Combining Equations (5-12), (5-13), (5-16), and (5-17) the 
following relationship is obtained 


dR 


R 


s dt 


A 

u A + 

SDL s 


i/2 


R 


= - 1 


(5-18) 


Separating variables and integrating Equation (5-18) , noting 


that T =0 at Rg = 1, yields 
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T = 


(1 - R ) 
s 

1/2 


A 


SDL 


In 


A 

X/ X 

m + 

R 

°SDL \\ ■ 

s 


m 


SDL 


. 1/2 


(5-19) 


where m and A are functions of pressure (A for the non-wet bulb 
SDL ® ® 

case is relatively constant with pressure for E = 10 and n = 2) and 

A^ is a given constant depending upon the injected droplet size. 


For the limit as A^/A^ approaches zero, integrating Equation 
(5-18) yields 


T 




SDL 


(5-20) 


which represents the size variation for a small droplet. 

Setting R == 0 in Equations (5-19) or (5-20) yields the 

dimensionless lifetime x for drops originating with a radius r *. 

So 

The fraction of these droplets in a size range R ’ - R " is then 

s s 

given by 


F(R^’ - (5-21) 

The response has been determined as a function of co and A, 
at a given pressure and injected liquid temperature. These quantities 


are related to values based upon the injected radius, cii , and A 
as follows 


A . = A R ^ , CO . = (u R 
J o s . ’ j os, 
' J 3 


(5-22) 


for any size group, j. The response of each size group, j, is then 
summed using Equation (5-9) to give the total response. A second 
weighted sum is required if more than one initial drop size is 
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i, 


f ■-* 

I 

\} 

1 = 


injected. In this case a reference size is more convenient and 



In the following, the reference size has been taken to be the largest 
droplet size injected, and the total response is calculated from 
Equation (5-11) . 




CHAPTER VI 


SPRAY COMBUSTION RESULTS 
6.1 Droplet Burning Time 

Employing Equation (S-9) , calculations of droplet lifetimes 

were completed for a range of conditions. Figure 26 shows the 

nondimensional droplet lifetime as a function of and pressure, 

for both the wet and non-wet bulb eases. At the large drop limit, 

the dimensionless lifetime is relatively independent of pressure. 

This is true for a second order reaction, where the burning rate per 

unit area is proportional to pressure. At the small drop limit for 

the non-wet bulb case, increases with increasing pressure, due to 

the reduction in m with pressure. For the wet bulb case T 

SDL L 

decreases with increasing pressure at the small drop limit, since 

increases with increasing pressure at this condition, 

SDTj 

Figures 27 and 28 show the variation of dimensionless radius 
with dimensionless time and A^/A^; at a pressure of 10 atmospheres 
for both non-wet and wet bulb conditions. These figures show that 
increases, decreases for a given pressure. This is due 
to the increased gasification rate caused by decomposition near the 
droplet surface. Results at 1 and 100 atmospheres are similar to 
Figures 27 and 28. 

Table 7 lists the dimensional droplet lifetimes as a function 
of pressure and initial drop radius. 


NON-WET BULB 



WET BULB 



Figure 27 


Nondimensional Drc 
of Nondimensional 
10 Atm. Pressure i 







T 


Figure 28 Nondimensional Droplet Burning Time as a Function 
of Nondimensional Droplet Radius at Approximatelv 
10 Atm. Pressure (Wet Bulb Conditions) 
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Table 7 


Mmensior.al Droplet Lifetimes as a Function of 
Initial Radius and Pressure 


p A 
(N/m2) 

t*(M) = 10 
so 

100 

1000 

(Non-Wet Bulb) 


(seconds) 


10^ 

0 .00158 

0.128 

4.79 

10^ 

0 .00135 

0.0489 

0.728 

10^ 

0.000507 

0.00734 

0.0788 

(Wet Bulb ) 


(seconds) 


10^ 

0.00141 

0.109 

3.61 

10^ 

0.000955 

0.0272 

0.361 

10^ 

0.000179 

0.0020 

0,0229 


The wet bulb droplets eonsistently have shorter lifetimes 
because of their higher mass burning rates. For all eases, at a given 
initial radius, the lifetime decreases with increasing pressure. 

This is caused by the fact that a given radius corresponds to a 
larger value of A at higher pressures, yielding a greater mass 
burning rate< At a constant pressure, the lifetime increases with 
drop size, which is to be expected. 

The response portion of the present analysis is not valid at 
frequencies having an oscillation period comparable to the lifetime 
of the droplet. This provides a minimum droplet size that can be 
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considered by the present analysis. The probable oscillation periods 
of instability range from 2 x 10“^ to 3 x 10“^ seconds/eycle 
(corresponding to the 500-30,000 Hz range indicated in Reference [1]). 
At the oscillation period limit of 2 x lO'^ seconds /cycle, the minimum 
drop radius is approximately 13 p at pressures of 1 and 10 atmospheres, 
for both wet and non-wet bulb conditions. At a pressure of 100 
atmospheres, the minimum drop radius is approximately 30 p at the 
non-wet bulb condition, and approximately 95 p for the wet bulb case. 

At the other limit of 3 x 10 seconds/cyde the minimum droplet 

radius is approximately 2 p at all pressures and both wet and non-wet 
bulb conditions. This indicates that at low frequencies, appreciable 
response may be generated by droplet lifetime characteristics, a 
mechanism that has not been considered in the present investigation. 

^ ^ Qpdsi~Steady Gas Hhase Spray Response 

If a monodisperse spray is injected into the combustion 
chamber, the total response, is the weighted sum of all drop 

sizes existing in the chamber, Pigure 29 illustrates the non-wet 
bulb total response, of a monodisperse spray as a function of frequency, 

for a pressure of 10 atmospheres. The values of A were chosen 

o 

so that the first drop size group (which has the largest mass 
fraetron) has an A for which the droplet response has previously 
been computed. The solid lines correspond to the entire spray while 
the dashed lines oorfespond to the situation when it is assumed that 
all drop sizes are approximated by a constant drop size equal to A . 

Per these two cases, the peaks occur at approximately the same 
frequency, while the magnitude of the peak is slightly larger when 
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Figure 29 Total Response as a Function of cOq for a 

Monodisperse Spray at Approximately 10 Atm^ 
Pressure (Non-^Wet Bulb Conditions) 
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the drop size is assumed constant. The reduction in the response 
peak of the spray is due to the reduced response of the smaller 
droplets. It should be noted that for the case when the drop size 
xs constant at « 1,234 x 10^®, transient gas phase effects begin 
to appear at approximately = 6.98 x 10^. This corresponds to 

u = lO"^. 
s 

Figure 30 is the wet bulb total response for a monodisperse 
spray as a function of oj^ at 10 atmospheres pressure. The same 
conclusion can be drawn as from the non-wet bulb case. The total 
responses at pressures of 1 and 100 atmospheres for both wet and 
non-wet bulb conditions give similar results. The shapes of the 
total response curves in Figures 29 and 30 are very similar to the 
shape of the Individual response curves for A^. In all cases, five 
drop size groups were used, in which approximately 60% of the total 
mass occurred in the first group. 

When a poly disperse spray is introduced into the combustion 
chamber, the spray is separated into discrete spray sets, each with 
its own A^. The smaller initial drop size spray sets are placed in 
the corresponding size groups of the largest spray set as shown in 
Figure 25. The responses are weighted and summed, and the total 
response is obtained. 

Figure 31 is a plot of a non— wet bulb polydisperse spray at a 

pressure of 10 atmospheres, consisting of three spray sets with A 's 

o 

of 1,234 X 10 , lx 10 , and 1.234 x lO*^ . The response from each 

size group is weighted and summed in order to obtain the total 
response. The mass fractions of the spray sets and injection rates 
are varied, as listed in Figure 31, The total responses are greater 
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Figure 31 Total Response as a Function of for a Poly disperse Spray at Approximately 10 Atm 
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than unity when the larger spray sets comprise about 30 % of the 
injected spray mass. The results indicate, however, that even a 
very small percentage of large drops in the injected spray can result 
in substantial response. Findings at other pressures for both wet 
and non-wet bulb conditions are similar; however, for the wet bulb 
ease, no peaks greater than unity are observed. 


CHAPTER VII 


SUMMARY AND CONCLUSIONS 

7.1 Siiiiimary 

The present study considered the steady state combustion and 
oscillatory response of a single monopropellant droplet subjected to 
Imposed pressure oscillations, followed by an estimation of the 
characteristics of an entire spray „ The specific objectives of the 
study are as follows; 

1. Develop a mathematical model for the burning of a liquid 
monopropellant droplet for both steady state and 
oscillatory pressure conditions. 

2. Theoretically investigate the steady state droplet 
burning for both wet and non-wet bulb conditions; and 
compare the results with available measurements. 

3. Investigate the theoretical burning response for wet and 
non-wet bulb conditions. 

3 . Develop a mathematical model for one-dimensional , quasi- 
steady spray combustion, and calculate the overall spray 
combustion response for both wet and non-wet bulb 
conditions . 

The present theoretical model is an extension of the one- 
dimensional strand model developed by Allison and Eaeth [11] . The 
major assumptions in the analysis are: a one— step. Irreversible 

gas phase reaction, conventional low pressure phase equilibrium at 
the liquid surface; constant liquid phase physical properties, and 
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variable gas phase physical properties. The present calculations were 
based on the reaction and the physical property characteristics of 
hydrazine. 

A limitation of the work involves neglecting forced convection, 
yielding a spherically symmetric flow field. It is believed that 
low droplet velocities represent a critical condition for instability 
[1] , and forced convection is of minor importance when the reaction 
zone is very near the surface [25, 26], a condition where the response 
for menopropellants is greatest. 

The theoretical approach employs the use of a perturbation 
analysis in which the amplitude of the imposed pressure oscillation 
is taken as the small perturbation parameter. 

The oscillatory (first-order) solution involves a transient 
liquid and both transient and quasi-steady gas phase effects . The 
results indicated that most of the transient gas phase effects lie 
outside the probable instability range at realistic combustion 
chamber pressures. Therefore, the spray response calculations were 
confined to conditions where a quasi-steady gas phase is valid. 

Finally, a one-dimensional quasi-steady spray combustion 
model is developed in which the spray is divided into size groups , and 
the response is weighted and summed over all groups , to give the 
total spray response. 

7.2 Conclusions 

This section is divided into three parts: steady state 

combustion, oscillatory combustion, and spray combustion. 
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7.2.1 Steady State Combustion . The major conclusions of the 
steady state portion of this study are as follows: 

1. The theoretical droplet surface temperaturos are identical 
to the measured and predicted strand values from Reference 
[ 11 ] . 

2. Steady state droplet mass burning rates can be divided 

into three regions: large, intermediate, and small drop 

regions . The large drop region can be approximated by a 

one-dimensional strand, while the small drop region is 

similar to a small drop evaporating in absence of reaction. 

Most drop sizes of technological importance fall within 

4 9 

the intermediate region (10 < A < 10 ) . 

3. At the non-wet bulb small drop limit, the burning 
rate decreases with increasing pressure because the 
higher pressures raise the surface temperature, while 
the ambient temperature is constant. The reduced 
temperature difference between the liquid surface and 
the ambient gas causes a reduction in the droplet 
gasification rate, 

4. The wet bulb small drop limit burning rate increases 
with pressure. This is caused by the increase in the 
ambient temperature with pressure, which compensates j >r 
the increase in surface temperature. 

5 . At all pressures , the wet bulb burning rate is higher 
than the non-wet bulb case. 'This is due to the fact that 
the wet bulb droplet has a higher ambient temperature, and 


from the fact that the heat conducted back into the wet 
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bulb droplet goes only into vaporizing the liquid, not 
in heating the droplet liquid. 

7.2.2 Oscillatory Combustion . The major conclusions of the 
oscillatory portion of this study are as follows: 

1. At the small drop limit, the response is negative and 
becomes more negative with increasing pressure. This 
results from the gasification rate decrease as the surface 
temperature increases for a constant ambient temperature. 
(The wet bulb response is less negative than the non-wet 
bulb case.) Therefore, pressures above the mean result 

in reduced gasification rates while pressures below the 
mean result in increased gasification rates - yielding a 
negative response, 

2. The large drop limit results give good agreement with the 
one-dimensional strand results of Reference [11] . 

3. At the non-wet bulb large droplet limit, response peaks 
occur near the characteristic frequency for the liquid 
and gas phase thermal waves. With increasing pressure, 
these liquid and gas phase peaks approach each other 
and increase in magnitude. 

4. At the wet bulb large drop limit, no response peaks 
occur from liquid transient effects, 

5. As the drop size decreases, the response also decreases. 

6. The transient gas phase response is not a major 
contributor to combustion instability at the frequencies 
and pressures where instability has a high probability 
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of occurring e Thus, a quasi— steady gas phase analysis 
can be assumed in most instances, 

7 . A response peak will occur only when the mean temperature 
gradient near the droplet surface has sufficient steepness, 
so that enough energy will flow from the surface towards 
the center of the drop, causing the perturbation 
temperature to lead the pressure. When the mean 
temperature gradient is not steep enough, (which occurs 
near the wet bulb case for properties used in the present 
calculations) the perturbation surface temperature will 
lag behind the pressure at all frequencies, 

7,2,3 Spray Combustion , The major conclusions of the spray 
combustion portion of the study are as follows: 

1. The total response of a polydisperse spray is less than 
the response of the largest droplet injected, 

2. Relatively small numbers of injected large drops greatly 
influence the spray by adding a substantial response, 

3. A spray response of order unity, due to the liquid phase 
thermal wave, has heen observed, indicating that this 

is a viable mechanism for producing combustion instability, 

4. At high pressures and low frequencies (on the order of 
500 Hz) the large burning rates of monopropellant dvroplets 
provide droplet lifetimes near the cycle period for 

drop sizes of technological importance. Under these 
conditions, droplet lifetime effectiS contribute to the 
response; thus, this phenomena should be investigated 


next , 
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